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Abstract 

We review the properties of the self-organized critical (SOC) forest-fire 
model. The paradigm of self-organized criticality refers to the tendency of 
certain large dissipative systems to drive themselves into a critical state inde- 
pendent of the initial conditions and without fine-tuning of the parameters. 
After an introduction, we define the rules of the model and discuss various 
large-scale structures which may appear in this system. The origin of the 
critical behavior is explained, critical exponents are introduced, and scaling 
relations between the exponents are derived. Results of computer simulations 
and analytical calculations are summarized. The existence of an upper crit- 
ical dimension and the universality of the critical behavior under changes of 
lattice symmetry or the introduction of immunity are discussed. A survey 
of interesting modifications of the forest-fire model is given. Finally, several 

other important SOC models are briefly described. 
PACS numbers: 05.40.-hj, 05.70.Jk, 05.70.Ln 
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I. INTRODUCTION 



Driven dissipative systems far from thermodynamic equilibrium show a rich variety of 
patterns. Since they receive a permanent input of energy, they can maintain states which 
are highly ordered or complex. Ubiquitous examples for complex structures are fractals [jl],0] 
as well as their temporal counterpart, 1/ f -noise [^Q. Fractals are self-similar structures 
that look the same on different scales of observation, since they have no intrinsic length 
scale. Their spatial correlation functions are power laws. Well-known examples are moun- 
tain landscapes and coastlines. l//-noise is the temporal equivalent of fractals. Its name 
indicates that the Fourier transform of the temporal correlation function is a power law of 
the form 1//° with a 1. Like fractals, l//-noise can be found in many natural systems, 
e.g. undersea ocean currents or stock market prices. 

In 1987, Bak, Tang, and Wiesenfeld introduced the sandpile model which evolves into a 
critical state irrespective of initial conditions and without fine tuning of parameters [^0. 
Such extended non-equilibrium systems are called self-organized critical (SOC) and exhibit 
power-law correlations in space and time. The concept of SOC (for a review see 0) has 
attracted much interest since it might explain part of the abundance of fractals and 1 / /-noise 
in nature and create a link between the two. Models for earthquakes the evolution of 
populations [|r^,^], the formation of clouds \T2\ and river networks |TB[, and many more have 
been introduced and investigated mainly by computer simulations, but with few exceptions, 
most of these SOC models still are barely understood. 

This review is mainly devoted to the forest-fire model (FFM) fl^ , which is a particularly 
simple example for open systems which shows nevertheless a variety of different structures, 
depending on the parameters. In the limit of a double separation of time scales, the FFM 
becomes SOC, with a power-law distribution of fires and forest clusters. 

The FFM is a stochastic cellular automaton which is defined on a (i-dimensional hyper- 
cubic lattice with L'^ sites. Each site is occupied by a tree, a burning tree, or it is empty. 
During one time step, the system is parallely updated according to the following rules 
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1. Burning tree — > empty site. 

2. Tree — > burning tree with probability 1 — (7 if at least one nearest neighbor is burning. 

3. Tree — > burning tree with probability / if no nearest neighbor is burning. 

4. Empty site — > tree with probability p. 

An important application of the FFM comes from its close relation to excitable media 
which comprise phenomena so different as spreading of diseases, oscillating chemical 



reactions, propagation of electrical activity in neurons or heart muscles, spiral galaxies, and 
many more (For a review on excitable systems see e.g. ||T6| , p!7| ). These systems essentially 
have three states which are called quiescent (= tree), excited (= burning tree), and refractory 
(= empty site). Excitation spreads from one place to its neighbors if they are quiescent. 
After excitation, a refractory site needs some time to recover its quiescent state. Since 
excitable systems are often deterministic, it is, however, not always possible to describe 
their behavior in the framework of the FFM, where some of the parameters are stochastic. 

Turning back to the FFM, we consider a system with arbitrary initial conditions. After a 
transition period the system approaches a steady state the properties of which depend only 
on the parameter values. If the system size L is large enough the properties of the steady 
state are also independent of the boundary conditions. Let pe, pt, and p/ be the mean 
density of empty sites, of trees, and of burning trees in the steady state. These densities are 
related by the equations 

Pe + Pt + Pf = l (1.1) 

and 

pf=Ppe- (1.2) 

The second relation says that the mean number of growing trees equals the mean number of 
burning trees in the steady state. Depending on the magnitude of the parameters, a variety 
of large scale structures arise. 



A. Fire fronts and spirals 



In the limit of slow tree growth p and without lightning and immunity, fire only spreads 
from burning trees to their neighbors but does not occur spontaneously. In this version, 
which was originally introduced by P. Bak, K. Chen, and C. Tang [0, the fire fronts 
become more and more regular and spiral-shaped with decreasing p . A snapshot of this 
state is shown in Fig. [I|. A more thorough analysis of the spatial and temporal correlations 
in the system revealed the existence of characteristic length and time scales proportional to 
both of which become increasingly distinct with decreasing J9 indicating increasing 



determinism in this limit. Spiral waves are very familiar in excitable media and have been 
observed in oscillating chemical reactions, the heart muscle, and the chemical signals of the 
amoeba dictyostelium discoideum ||T6| , p!7|| . It is also known that certain epidemics occur 
periodically [0. They occur already in simple deterministic versions of the FFM 



B. Percolation due to immunity 

When the immunity is nonzero, the fire fronts present for (7 = become more and more 
fuzzy with increasing immunity g, and the forest becomes denser [^. At a critical value 
g = gc{p) the fire density becomes zero and the forest density becomes one. Figs. |^ and |^ 
show two snapshots of the system for values of g far below gc{p) and near gc{p)- At gc{p), 
the fire just percolates through the system. Since sites are not permanently immune, this 
kind of percolation is different from usual site percolation and is in the same universality 
class as directed percolation in d + 1 dimensions (the preferred direction corresponds to the 
time) [^. The critical immunity gc{p) increases with increasing p, since the fire then can 
return sooner to sites where it already has been. Above gc{p) the steady state of the system 
is a completely dense forest. A similar model, using the language of spreading diseases has 



been studied independently in |23 . 
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C. Self-organized critical (SOC) behavior 



The SOC behavior occurs when the hghtning probabihty is nonzero. For simphcity, we 
set the immunity to zero, but we will show later that the SOC state persists for g > 0. The 
ratio p// is a measure for the number of trees growing between two lightning strokes and 
therefore for the mean number of trees destroyed per lightning stroke. In the limit 

(1-3) 

there exist consequently large forest clusters and correlations over large distances. The 
model is SOC when tree growth is so slow that fire burns down even large clusters before 
new trees grow at their edge. This condition, which assures invariance under a change of 
the length scale, reads 

p<T-i(w), (1.4) 

where T(smax) is the time the fire needs to burn down a large forest cluster. It diverges in 
the limit f /p — > according to T(smax) oc {f /p)"^', with an appropriate exponent v' . 

In this case, the dynamics of the system depend only on the ratio //p, but not on / 
and p separately. When / and p are both decreased by the same factor, the overall time 
scale of the system is also changed by this factor, but not the number of trees that grow 
between two lightnings and therefore not the size distribution of forest clusters and of fires. 
In the simulations, the condition ( |1.4| ) is most easily realized by burning forest clusters 



instantaneously, i.e. during one time step. This extreme limit of the SOC forest-fire model 
has been invented independently by C. Henley 



The inequalities ( |1.3[ ) and ( |1.4| ) represent a double separation of time scales 

uipr' ^p-' (1.5) 

which is the condition for SOC behavior in the FFM. The time in which a forest cluster burns 
down is much shorter than the time in which a tree grows, which again is much shorter than 
the time between two lightning occurrences. Separation of time scales is quite frequent in 



nature, while the tuning of parameters to a certain finite value only takes place accidentally. 
Thus, the FFM is critical over a wide range of parameter values. A snapshot of the critical 
state is shown in Fig. ^. 

In the following, we will mainly discuss the properties of the SOC FFM, and we will 
compare this model to other SOC systems, as there are e.g. the sandpile model, the earth- 



quake model, etc. The outline of the remainder of this paper is as follows: In Sec. p A|, the 
critical exponents of the forest-fire model are defined, and scaling relations between them 
are derived. Sec. [II B| presents the results of computer simulations and discusses the issue 
of universality of the critical behavior. Sec. |II C| deals with the analytical results including 



renormalization group approaches. In Sec. Ill, several modifications of the forest-fire model 



are discussed. Some of the most prominent other SOC systems are presented in Sec. |V. 
Sec. contains concluding remarks. 



II. THE SELF-ORGANIZED CRITICAL FOREST-FIRE MODEL 
A. Critical exponents and scaling relations 

In this section, we consider some principal properties of the SOC FFM which will lead us 
to the definition of critical exponents and the derivation of scaling relations between them. 
In the following, we restrict ourselves to the case g = 0. 

The mean number s of trees that are destroyed by a lightning stroke can be calculated 
as follows. During one time step, there are fptL'^ lightning strokes in the system and ppeL'^ 
growing trees. In the steady state, the number of growing trees equals the number of burning 
trees, and therefore the mean number of trees destroyed by a lightning stroke is 

, = m^El^, (2.1) 

fpt f pt 

In the last step, the fire density which is very small due to time scale separation was ne- 
glected. For small values of f /p, the forest density pt approaches a constant value. If this 
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constant value is less than 1, the second factor on the right-hand side of Eq. (|2.1| ) is also 
constant for small //p, and Eq. (|2.1|) then represents a power law 



s oc if/p)-\ (2.2) 
In d > 2 dimensions, the critical forest density 



Pt = 1/m Pi 

f/p^o 



in fact, must be less than 1, as the following consideration indicates: If the critical forest 
density were = 1 in d > 2 dimensions, pt would be very close to 1 for small values of f/p. 
Then the largest forest cluster would contain a non-vanishing percentage of all trees in the 
system, and the average number s of trees burned by a lightning stroke would diverge in the 
limit L — > oo with fixed f/p, in contradiction to Eq. ( |2.1| ). In one dimension, there is no 
infinitely large forest cluster in the system as long as p^ < 1 , and therefore the critical forest 
density is p^ = 1. Nevertheless Eq. ( p.2| ) holds also in 1 dimension since the forest density 
approaches its critical value only logarithmically slowly, as will be shown in Subsection [11 C| . 
Eq. ( |2.2| ) indicates a critical point in the limit f/p 0. Close to this critical point, i.e. if 
f p, there is scaling over many orders of magnitude. 

Let n{s) be the mean number of forest clusters per unit volume consisting of s trees. 
Then the mean forest density is 



oo 



Pt = J2sn{s), (2.3) 
1 

and the mean number of trees destroyed by a lightning stroke is 

oo 

s = Y.sMs)/Pt. (2.4) 

1 

Since lim//p^o Pt is finite and s diverges as {f /p)~^, these equations imply that n{s) decreases 
at least like but not faster than s~^. As long as the system is not exactly at the 
critical point f/p = 0, i.e. for non- vanishing f/p, there must be a cutoff in the cluster size 



distribution for very large forest clusters. We conclude that |jT4| 



n{s) OC S ^C(s/Smax) (2.5) 

with 2 < r < 3 and 

Smax OC (//p)"^ CX s\ (2.6) 

The cutoff function C{x) is essentially constant for x < 1 and decreases to zero for large x. 
Eqs. ( |2.4| ) - ( |2.6|) yield s oc sj^"^' which leads to the scaling relation 

A = l/(3-r). (2.7) 



In the case r = 2, the right-hand side of Eq. (|2.5|) acquires a factor l/ln(sr 
the forest density given by Eq. (|2.3|) must not diverge in the limit f/p 0. The mean 
number of forest clusters per unit volume ^i^), therefore, decreases to zero for f/p — > 0, 
and consequently the forest density approaches the value 1. This situation occurs in one 
dimension. 

We also introduce the cluster radius R{s) (radius of gyration) which is the mean distance 
of the trees in a cluster from their center of mass. It is related to the cluster size s by 

s OC i?(s)^, 

with the fractal dimension fi. 

The correlation length ^ is defined by 

We conclude 

^ OC {f/p)-" with V = A//i. (2.8) 

In percolation theory, the hypers caling relation 

d = /i(r-l) (2.9) 

is satisfied, but it is not satisfied in the SOC FFM in = 2, as first stated in where 
also an interpretation of this relation is given: If Eq. ( p.9|) is satisfied, every box of Z'^ ^ 1 



sites contains a spanning piece of a large cluster when the system is at the critical point. In 
the FFM, there are at least in d = 2 many regions which contain no large forest cluster (see 
Fig. H), and consequently d < fi{T — 1). 

The mean forest density pt approaches its critical value = lim^/p^o pt via a power law 

Pi-Pt'xif/py/'. (2.10) 

One can also introduce several dynamical quantities and corresponding exponents char- 
acterizing the temporal behavior of the fire (see [^). Here, we mention only the time 
T{saia.x) oc {f /p)'"' which it takes to burn down a cluster of size Smax- 



B. Simulation results 

In this subsection, the simulation results for the critical exponents of the FFM are 
presented and discussed. 



1. Critical Exponents 

In one dimension, the critical exponents were determined not only by computer simula- 
tions, but also analytically |2^ (See also Subsection |II C 2|) . In higher dimensions one has 



to resort to computer simulations. The model is simulated most effectively with a method 
first proposed in pSl, which iterates the following rules: 



1. Choose an arbitrary site in the system. If it is not occupied by a tree, proceed with 
rule 2. If it is occupied by a tree, then ignite the tree and burn down the forest cluster 
to which the tree belongs. While burning the trees, evaluate the properties of the 
cluster as size, radius, etc. Proceed with rule 2. 

2. Choose pI f arbitrary sites in the system and grow a tree at all chosen empty sites. 
Proceed with rule 1. 
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By these rules, time scale separation is perfectly realized, and Eq. is satisfied. In order 
to assure that the system is in the steady state, a sufficiently large number of time steps 
have to be discarded in the beginning of each simulation. In the following, we discuss the 
simulation results obtained for the critical exponents and the critical forest density. 

The exponent r is determined by measuring the cluster distribution n{s) or the fire dis- 



tribution sn{s), which gives r — 1. In c? = 2, r ~ 2.15, as found in ||28|j29|j25|j26|| . With 



increasing dimension, r increases, too |2^j2^, and assumes its mean-field value 2.5 in di- 
mensions larger than 6 (for simulation results see Fig. ^, and for a complete list of the 
values see Tab. 1). 

Reliable results for the exponent A can be achieved by considering the normalized in- 
tegrated distribution function P{s) = ds's'n{s'), which was first introduced in ||2^. By 
collapsing different curves of P(s) for different values of f /p , A is found to be 1.16 p5| , p6| . 

The critical forest density is ~ 0.41 in d = 2 p8| , |29|p5| , |26[| , and it decreases with 
increasing dimension p9|j26|] . Tab. | lists the results obtained in where the system size 
was larger than in [^|. The exponent l/(5is~0.5in(i = 2 [^,^,^. 

The fractal dimension /i of the forest clusters is obtained from the slope of the cluster 
radius R{s) (see Fig. In 2 dimensions fi ~ 1.96 P3| , |2B[] , which is smaller than 2, in contrast 
to earlier assumptions ||l4| , p8| ). The values in higher dimensions are given in Tab. I PI and 
seem to approach the mean-field value 4, which is supposedly exact above 6 dimensions. The 
hyperscaling relation Eq. 



is definitely violated in d = 2 (see for an interpretation), 
but cannot be ruled out from simulation results in higher dimensions. 

The correlation length ^ is dominated by large clusters and consequently large radii. 
Therefore, the exponent u could be directly determined only in d = 2 dimensions, with the 
result u ~ 0.58 p8| , |25|p6| . The temporal analogon to the correlation length is T(sinax) oc 
if/p)-"', with z/' ~ 0.61 in d = 2 [gg. 

All these simulation results for the critical exponents suggest that the SOC forest-fire 
model has an upper critical dimension dc = 6, above which the critical exponents are iden- 
tical with those of mean-field-theory, which again is identical to the mean-field-theory of 
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percolation . The strongest evidence for this behavior comes from the exponent r, which 



approaches the percolation value Tpcrc = 5/2 for d ^ 6 and is indistinguishable from 5/2 
in all simulated dimensions d > Q. But also in the other exponents the difference between 
forest-fire and percolation values seems to vanish with increasing dimension. 

2. Universality 

In equilibrium systems, the critical behavior usually depends only on properties as di- 
mension and conservation laws, but not on microscopic details. It can therefore be expected 
that also the critical behavior of the SOC FFM is universal under certain changes of the 



model rules. This assumption was checked in [^, where the 2D simulations were carried 
out for a triangular lattice {pi ~ 34%) and for a square lattice with next-nearest-neighbor 
interaction (p^ ^ 28%). The critical exponents of these variations of the model were found 
to be the same. See, for instance. Fig. |^, where the correlation length for the three models 
is plotted. The critical densities are smaller than for the square lattice, since the fire has 
more possible paths to spread due to the larger number of neighbors. 

3. Immunity 

Another modification of the model rules is obtained by including a non zero immunity 
g [|30|. In contrast to the rule given in Sec. |, the simulations were performed with immune 



bonds instead of immune sites. Not all trees that are neighbors of a burning tree catch fire, 
and consequently the fire does no longer burn forest clusters but clusters of trees that are 
connected by non immune bonds. A small value of immunity corresponds essentially to a 
change in the lattice symmetry, since the effective coordination number is decreased. As 
we have just seen, a change of the lattice symmetry does not affect the critical exponents, 
which is again confirmed by the simulations. With increasing immunity, the forest density 
increases until the critical forest density becomes = 1 at a critical immunity Qc- Here, 
a new scenario occurs: The forest is completely dense in the limit f/p — > 0, and clusters 
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that are destroyed by fire are percolation clusters of bond percolation. Consequently, the 
exponents Tc = T{g = g^) and are given by percolation theory, and the threshold is 
Qc = 1/2, which is 1 minus the bond percolation threshold. Since the critical forest density 
is 1 at (7 = Qc, scaling relation Eq. ( |2.7| ) now reads 

Ac = 7c/(3-r,) (2.11) 

with a new exponent 7c. 

For finite f /p, the mean forest density is no longer 1, and one can determine also the 
other critical exponents Ac, 5c, i^c, and 7c. The scaling relations Eq. (|2.11]) and Eq. ( p.8|) 
are confirmed by the simulations. Using Eqs. ( |2.1j ) and ( p.lO|) , one obtains a new scaling 



relation s oc {f /p)^^" with 7c = 1 — l/5c, which is also confirmed by the simulations. 

When the immunity is just below its critical value Qc-, a crossover from percolation-like to 
SOC behavior is observed. On length scales smaller than the percolation correlation length 
■Cperc OC ((7c ^ fi')^''"^ & systcm close to the percolation threshold cannot be distinguished 
from a system exactly at the percolation threshold. As long as //p is so large that fires do 
not spread further than ^perc, the exponents are identical to those aX g = Qc- When f /p 
becomes very small, there are fires which spread further than the percolation correlation 
length. These fires are stopped by empty sites that were created by earlier fires. This is the 
same mechanism as for small g: fires that would spread indefinitely if there were no empty 
sites are stopped by empty sites. We conclude that these large fires lead to the critical 
exponents A, z/, and 5 that have been observed for = 0. We make the following scaling 
ansatz for the correlation length: 

e = (//pr^F (^) . (2.12) 

It is plausible that the crossover from percolation-like to SOC behavior takes place when 
f /p becomes so small that the correlation length exceeds the percolation correlation length, 
which suggests that the crossover exponent is given by = z/c/ Vperc- The scaling function 
F{x) is constant for small x and is oc x^'^''^"'^/'^ for large x. Analogous scaling laws hold for 
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Smax and pl — pt- Fig. ^ shows the simulation results for the scaling function of the correlation 



since all curves coincide. The dashed line represents -F(O) as obtained from the simulations 
at Qc- This crossover is similar to crossover phenomena at equilibrium phase transitions, and 
it should be observed also in higher dimensions. Although the crossover in r and /i vanishes 
in dimensions larger than 6, there is still a crossover in A and due to the modified scaling 
relation at Qc- In = 1, the critical immunity is Qc = 0, and no crossover can take place. 



In this subsection, we review several analytical approaches to the SOC forest-fire model. 
First, we present the mean-field theory (MFT) which leads to the same exponents as the 
MFT of percolation. Then, we give several exact results for the one-dimensional model. 
Finally, we discuss attempts to establish a renormalization group for the forest-fire model. 



One of the simplest possible analytical approaches to a model is a MFT, which neglects 
all correlations in the system and describes it entirely in terms of densities. The neglect 
of correlations can be modelled in simulations by constructing a random-neighbor model, 
where the nearest-neighbor connections are chosen randomly at each time step p9|. The 
probability that a randomly chosen neighbor of a given burning tree is occupied by a green 
tree is obviously given by the tree density. In the stationary state, the mean-field equations 



length F{x) for different values of Qc — g. The scaling ansatz Eq. ( p.l2| ) is well confirmed 



C. Analytical results 



1. Mean- field theory 



are the following |2|,|3T]j3|: 



Pf = PPe , 



(2.13) 




(2.14) 



Pe + Pt + Pf = I- 



(2.15) 
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The first and last equation are exact (see Eg. ( |1.1|) and Eq.( |1.2| )), since they involve no 



nearest-neighbor interactions. The second equation includes a mean-field approximation, 
since the probability that one or more neighbors of a given site are burning is given by 
(1 — (1 — pfY'^), without taking into account any correlations. 

From these equations, we find an implicit equation for the fire density alone 

P/ = (1 - Pf (1 + Up)) (i - (1 - /) (1 - Pf)'') ■ (2.16) 

In the SOC limit, p, /, and pf are very small compared to 1, and we find to leading order 

Pf = p{l - l/2d) + f/2d + 0{p',pf, ■■■), (2.17) 
Pe = l- l/2d + f/2dp + 0{p, /,•■■), (2.18) 
pt = l/2d-f/2dp + 0{pj,---). (2.19) 

The tree density approaches its critical value linearly in f /p, leading to 6 = 1. The critical 
tree density is Pt{f/p = 0) = l/2d, which means that a burning tree ignites on an average 
just one other tree. This situation is identical to MFT of percolation or, equivalently, 
percolation on a Cayley tree, where percolation proceeds to each available neighbor with 
the same probability and to one neighbor on an average. The cluster size distribution and 
the fractal dimension of clusters for this problem are well known We therefore obtain 



T = 2.5 and /i = 4. Scaling relations Eqs. ( p.7|) and ( p.8|) then give A = 2 and u = 0.5. 
Percolation has the upper critical dimension 6, and one might expect that the SOC forest- 
fire model has the same upper critical dimension, as conjectured in [^,^,^] and supported 
by the simulation results (see subsection PB| ). 

When the limit p ^ with / = is considered instead of the limit / ^ with p <C f/p, 
the MFT of the forest-fire model also gives a critical tree density l/2d, leading again to the 
same exponents as MFT of percolation. In MFT the structural information is lost and thus 
MFT cannot see the qualitative difference between the SOC state and the quasideterministic 
state with spiral-shaped fire fronts. 
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2. Exact results in one dimension 



In one dimension, exact results have been obtained in |27,31,3i], thus proving analytically 
that non conservative models can indeed show SOC. Here, we give an intuitive derivation 
of the results: Consider a string of k p/f sites. This string is too short for two trees to 
grow during the same time step. Lightning does not strike this string before all of its trees 
are grown. Since we are always interested in the limit f /p ^ 0, the following considerations 
remain valid even for strings of a very large size. Starting with a completely empty string, 
it passes through a cycle which is illustrated in Fig. ^. During one time step, a tree grows 
with probability p on any site. After some time, the string is completely occupied by trees. 
Then the forest in the neighborhood of the string will also be quite dense. The trees on 
the string are part of a forest cluster which is much larger than k. Eventually that cluster 
becomes so large that it is struck by lightning with a non-vanishing probability. Then the 
forest cluster burns down, and the string again becomes completely empty. 

This consideration allows us to write down rate equations for the states of the string. 
In the steady state, each configuration of trees is generated as often as it is destroyed. Let 
Pjt(m) be the probability that the string is occupied by m trees. Each configuration which 
contains the same number of trees has the same probability. A configuration of m trees is 
destroyed when a tree grows at one of the empty sites, and is generated when a tree grows 
in a state consisting of m — 1 trees. The completely empty state is generated when a dense 
forest burns down. Since all trees on our string burn down simultaneously, this happens 
each time when a given site of the string is set on fire. This in turn happens as often as a 
new tree grows at this given site, i.e. with probability p (1 ~ pt) per time step. We therefore 
have the following equations, 

pkPk{0)=p{l-pt), (2.20) 
p{k — m)Pk{m) = p(k — m + l)Pk{m — 1) for m 7^ 0, A; . (2.21) 

We conclude 
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Pf,{m) = {1 - pt)/{k -m) for m< k, (2.22) 

k—l 

P^{k) = l-{l-pt)J2l/{k-m) 

m=0 
k 

= l-{l-pt)Y.l/m. (2.23) 

m=l 

These last two equations contain a wealth of information: cluster size distribution, hole 
distribution, growth velocity etc. 

A forest cluster of size s is a configuration of s neighboring trees with an empty site at 
each end. The size distribution of forest clusters is consequently given by 

Ps+2is) 1-pt 2 .2 24) 

(-+2^) (s + i)(s + 2) P''' ■ ^ ' 

This is a power law with the critical exponent r = 2. The size distribution of fires is 
oc sn{s) oc s~^. Fig. ^ shows the numerical result for the fire distribution sn{s). It agrees 
perfectly with Eq. ( p.24| ) in the region s < Smax- 

Sraax, iutroduccd in Eqs. ( |2.5| ) and ( p.6| ) is the characteristic size (in one dimension length) 
where the power law n{s) oc breaks down. We calculate Smax from the condition that a 
string of size k < Smax is not struck by lightning until all trees are grown. When a string of 
size k is completely empty at time t = 0, it will be occupied by k trees after 



T{k) = (l/p) ^ 1/m ~ ln(A;)/p 

m=l 

time-steps on an average. The mean number of trees after t time-steps is 

m(t) = k[l — exp(— pt)]. 
The probability that lightning strikes a string of size k before all trees are grown is 

T(fc) 

/ m{t) ^ U/v)KHk) -i)^ Ulv)kHk). 

t=i 

We conclude 

Smax In ( Smax) OC p// for large p/ f , 
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leading to A = 1. 

Next we determine the relation between the mean forest density pt and the parameter 
f /p. The mean forest density is given by 



s=l 



;i-a)E 



+ !)(. + 2) 

;i -pt)ln(Smax)- 



Thus 

pt 



In(smax) - ln(p//) for large p/f. 



I-Pt 

The forest density approaches the value 1 at the critical point. This is not surprising since 
no infinitely large cluster exists in a one-dimensional system as long as the forest is not 
completely dense. 

The exponents characterizing the size distribution of forest clusters and fires remain the 
same when the fire is allowed to jump over holes of several empty sites [|3^, proving the 



universality of the critical exponents. A change in the range of the interaction corresponds 
to a change in lattice symmetry in higher dimensions, where universality also has been found 
(see subsection [11 Bp . 

The values of the exponents change, however, when the rule for tree growth is modified. 
When tree growth becomes more deterministic, the critical behavior eventually breaks down, 
and the system becomes synchronized [^. To illustrate this, let us look at the case of 



completely deterministic tree growth, where an empty site turns to a tree exactly T time 
steps after the site has become empty. Neighboring sites both of which happen to be occupied 
by a tree burn down during the same fire. They consequently turn to trees simultaneously 
and burn down together simultaneously for all future times, i.e. they are synchronized. 
After some time, the system will consist of synchronized blocks which are so large that they 
are struck by lightning before the neighboring blocks become occupied by trees. 
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As an intermediate case, let us next consider a tree growth rule where the life time 
distribution of empty sites is a constant 1/T over a time T. An initially empty string 
of k sites then grows the same number of trees during each time step, irrespective of the 
number of occupied sites. The string consequently becomes occupied after a finite time. The 
above-mentioned deterministic model is just a coarse-grained version of this model, and we 
therefore expect that the system is dominated by large fires. Due to stochastic tree growth, 
however, there exist also small forest clusters. At small scales k, we find that Pk{fn) has 
the same value for all m < k, leading to n{s) oc s"^. The simulation results Fig. |ll] show 
this power law for small cluster sizes, and the expected peak at large clusters, indicating 
synchronization. When the distribution of life times of empty sites has a power-law tail, the 
system shows SOC behavior with an exponent which depends on the exponent characterizing 



the tail 36 



3. Renormalization group approach 

In equilibrium critical phenomena, the invariance of the system under a change of length 
and time scale together with a partial trace operation leads to a renormalization group 
(RG) approach. A RG not only yields values (possibly approximate, if one has to use 
perturbation theory) for the critical exponents in low dimensions, but also information on 
the upper critical dimension and on universality classes. 

The separation of time scales and the avalanche structure of the forest-fire model and 
other SOC models make it difficult to find an appropriate RG formalism for these systems. 
Two approaches which have been taken so far ||37|j38[l do not yet discriminate between the two 
limits limp^olim/^o and lim/^o with p <^ {f /pY' ■ The real-space RG in [0, which neglects 
correlations, yields the correct exponents in one dimension, and good approximations in two 



dimensions. In I^Sj, the FFM is mapped on a field theory, but the subsequent calculations 
contain severe approximations. In a third approach |3^, the forest-fire model is also mapped 
on a field theory. However, the renormalization in the limit of double time scale separation 
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remains a problem. 



III. MODIFICATIONS OF THE FOREST-FIRE MODEL 

The forest-fire model can be modified in many ways, and some of these modifications 
show very interesting behavior. A version with conserved tree density is suggested in |31[ 



and studied in |Q . After each fire, all trees which have been burnt are regrown at randomly 



chosen empty sites. For tree densities smaller than p^, this model has only finite fires and 
can be mapped on the SOC forest-fire model by finding the corresponding value of f /p. For 
densities between pi ~ 0.41 and a second critical density pf' ~ 0.435, the system shows 
critical behavior with a power-law size distribution of fires and of tree clusters. Both the 
cutoff in cluster size and the correlation length diverge with some power of the system size 
(see Fig. [T2|). The exponents depend on the density. Thus one has the interesting situation 
that criticality is not confined to a point but exists in a finite intervall. For even higher 
average density, the system splits into several subphases with different densities (Fig. p!3[ ). 
The subphase with highest density contains an infinite cluster. The number of subphases 
which can be sustained depends on the system size, their shape on the boundary conditions. 
Another modification of this model, where each tree is grown immediately after it has been 
burnt, shows both critical behavior and spiral-shaped fire fronts pl| . 

A deterministic version of the forest-fire model with continuous tree growth has been 



suggested in |^2|. The variable "tree height" is increased globally and very slowly. Trees 
higher than a certain threshold catch fire and ignite all neighbors above a second, smaller, 
threshold. The height of the tree after the fire is a function of the height before the fire. |^ 



studies a generalized version and shows that this model is not critical, but shows periodic 
behavior or finite fires, depending on the value of a parameter. Another continuous model, 
which additionally includes energy diffusion, shows small or large fires, depending on the 
parameters, and the two regimes are separated by a critical point. ^ . 
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IV. A SELECTION OF OTHER SOC MODELS 



Many systems have been introduced as examples for SOC. For part of these systems, the 
evidence is based on numerical data for relatively small system sizes with little analytical 
understanding of the origin of the critical behavior. In this section, we discuss three of 
the most thoroughly studied SOC models and compare them with the forest-fire model. 
Although the term SOC has been used in many different contexts, we restrict ourselves to 
systems with slow driving (energy input) and dissipation events which are instantaneous 
on the time scale of driving. The size distribution of dissipation events obeys a power law. 
Further examples for SOC include the dynamics of magnetic domain patterns |^ , ^ , cloud 
formation [|l^], evolution of populations |[lO|JTT| , erosion [jl3|, fracture fragmentation 



the Bean critical state in type II superconductors |^9[, depinning transitions of interfaces, 
charge density waves or superconducting flux lattices |^0|5 ^^"^ many more. Although not 
treated in this section, these examples underline the importance of SOC as a crossdisciplinary 
subject with applications in physics, biology, chemistry and geology. 



A. Sandpile model 

The sandpile model is the prototype for SOC Sand grains are dropped at random 

on the sites of a lattice. When the number of grains on a lattice site ("height model") 
or the difference in the grain number between neighbors ("slope model") exceeds a certain 
threshold, the site topples and redistributes its grains amongst the nearest neighbors. When 
a neighbor thereby is lifted above the threshold, it topples too, and the avalanche continues 
until all sites are below or at the threshold. Then a new grain of sand is added at random. In 
the stationary state, the mean number of added grains must equal the mean number of grains 
leaving the system at the edge, and there must exist avalanches spreading over a distance 
of the order of the system size. The sandpile therefore is in a critical state where sand 
grains are redistributed on all scales. For the height model, several exact results have been 
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obtained The critical state is robust with respect to a variety of changes, but is usually 
destroyed when the local conservation law is violated. Real sandpiles behave differently, 
mainly because large avalanches have large inertia and are not immediately stopped when 
the slope becomes small. 

In contrast to the forest-fire model, the sandpile model has only two separated time scales. 
The condition that avalanches relax before new grains of sand are added, corresponds to 
the condition that a fire is extinguished before new trees grow. In the sandpile model, 
large avalanches occur due to the local conservation of sand grains, in the forest-fire model 



large fires occur in the limit f /p 0- It has been pointed out in |^ that the sandpile 
model can be simulated exactly at the critical point, while the forest-fire model is only 
close to its critical point and has therefore a relevant parameter. However, in any realistic 
physical situation, the driving rate for the sandpile model can never be exactly equal to 
zero, resulting also in a relevant parameter and an upper cutoff in avalanche size. Above 
that size, avalanches overlap. 



B. Earthquake model 

The earthquake model p| is a continuous, non conservative variant of the sandpile model. 
It can be derived from a 2D spring-block model for the motion of tectonic plates (Burridge- 
Knopoff model). The force on all sites ("blocks") is increased uniformly and very slowly. 
When the force at a given site exceeds a threshold, the block moves to its equilibrium 
position, and the force on the corresponding site is reset to zero. The force on the 4 nearest 
neighbors is increased by a fixed percentage of the released force. This percentage depends 
on the ratio of the spring constants in the model, and is usually smaller than 1/4, i.e. the 
model has no conservation law. The system is not driven as long as there are forces above 
the threshold and an "earthquake" is going on. The size distribution of earthquakes is 
found to obey a power-law, in agreement with the famous Gutenberg-Richter law for the 
size distribution of real earthquakes . The exponent depends continuously on the degree 
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of conservation, and the critical behavior seems to persist even in the hmit where almost 
no force is transmitted to the neighbors. The size of the largest avalanche diverges with 
diverging system size, but so slowly that it covers only a vanishing part of the system 



5^ , |53[] . When the global driving is replaced by local driving, when inhomogeneities such 
as defects are present, or when the open boundary conditions are replaced by boundary 
conditions, the critical state is destroyed The authors of [|3 show that the critical 

behavior is related to a partial synchronization of neighboring sites. However, the origin 
of the critical behavior of the earthquake model, and in particular the dependence of the 
critical exponents on the degree of conservation, is relatively poorly understood. A similar 
behavior has been found so far only in one other model, namely the above-mentioned version 



of the forest-fire model with tree conservation ||4CI|| . 

In this context, another earthquake model that is also derived form the Burridge-Kopoff 
model, is worth mentioning ||5^. Like the above model, it is deterministic and continuous, 
but it includes additionally effects of inertia. It shows a power-law size distribution over a 
certain range of earthquake sizes, but the very large earthquakes occur more frequently than 
expected from a power-law distribution, and in periodic time intervals. From the evaluation 
of earthquake data, it is not obvious, which of the two models is more appropriate. 

C. Evolution model 



A particularly simple SOC model is the so-called "evolution model" [|Il|]. Each site in a 
one-dimensional chain is assigned a random number between and 1 ("fitness"). The prob- 
ability for a "mutation" depends exponentially on the fitness, and in the zero-temperature 
limit always the site with smallest fitness is mutated. This site is assigned a new random 
number between and 1. Since a mutation at one site affects the fitness of the neighbors, 
the two nearest neighbors are also assigned new random numbers. After some time, only 
a vanishing percentage of all sites have a fitness below a threshold ~ 0.67. A mutation of 
a site at the threshold releases an avalanche of mutations which is stopped when no site is 
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below the threshold any more. The size distribution of these avalanches obeys a power law. 
A review of the properties of this and related models and of some exact results is given in 

Reference mentions also a less spectacular application of this model, namely the 
depinning of interfaces, which usually occurs when a driving force exceeds a certain threshold. 
When not the driving force, but the interface velocity is chosen to be the parameter, the 
depinning transition becomes SOC in the limit of zero velocity |5^^8i, which is again a 
separation of time scales. 



V. CONCLUSIONS 

The main part of this review was devoted to the forest-fire model as a simple example 
for driven, dissipative systems with many degrees of freedom. Depending on the value of 
the parameters, the model shows fundamentally different structures. One finds a quasi- 
deterministic state with spiral waves, percolation- like behavior and, in particular, a self- 
organized critical state in the limit of double time scale separation. The properties of the 
SOC state have been analyzed numerically and analytically. The critical behavior has been 
shown to be very robust with respect to a variety of changes of the rules of the model. 

The notion of SOC has originally been introduced as possible explanation for the ubiquity 
of fractal structures and l//-noise in nature. Although this is a fascinating hypothesis, SOC 
seems to account only for the existence of part of these phenomena. Other extensively studied 
mechanisms producing fractals are diffusion- limited aggregation (DLA), kinetic roughening, 
and chaos and turbulence. 

All SOC models presented in this paper have slow driving and avalanche-like dynamics. 
However, not all systems that show these two features are SOC. Besides a power-law size 
distribution of avalanches, such systems might also have many small avalanches which release 
only little energy, or only large avalanches which release a finite part of the system's energy, 
or some combination of both. SOC systems are naturally at the critical point, due to 
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a conservation law (sandpile model), a second time scale separation (forest- fire model), 
a competition between open boundary conditions and the tendency of neighboring sites 
to synchronize (earthquake model), or due to the slow driving alone (evolution model). 
However, the critical behavior often breaks down when details of the model rules are changed, 
and is replaced by some other scenario. As examples we mentioned the forest-fire model with 
deterministic tree growth and the earthquake model with modified boundary conditions. 

The forest-fire model is closely related to excitable media. So far, spiral waves, target 
patterns and chaotic behavior have been found and investigated in excitable media. By 
investigating the appropriate parameter region one should also be able to find percolation- 
hke behavior and the self-organized critical state, which exist in the FFM. The SOC state 
should occur whenever there is spontaneous excitation that spreads very fast compared to 
the recovery time (e.g. fatal diseases, that occur seldom, but spread rapidly), provided that 
the probability distribution for the recovery time is not too narrow. 

After many years of studying formation of structure in non-equilibrium systems, one has 
got a glimpse at the mechanisms producing the overwhelming variety and complexity of 
structures surrounding us in nature, and there are certainly still many exciting phenomena 
to be discovered. 
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FIGURES 




FIG. 1. Snapshot of the Bak et al. forest-fire model in the steady state for p = 0.005 and 
L = 800. Trees are grey, burning trees are black, and empty sites are white. 
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FIG. 2. Snapshot of the forest-fire model far below the critical immunity for g 
and L = 200. Trees are grey, empty sites are white, and burning trees are black. 



= 0.2, p = 0.05, 




FIG. 3. Snapshot of the forest-fire model near the critical immunity for g = 0.48, p = 0.05, and 
L = 200. Trees are grey, empty sites are white, and burning trees are black. 
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FIG. 4. Snapshot of the SOC state in 2 dimensions. Trees are black, empty sites are white. 
The parameters are L = 1024 and f /p = 1/500. 
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FIG. 5. Normalized cluster size distribution n(s)/n(l) for d = 2 to 6 dimensions. The values 
of //p are 1/32000, 1/2000, 1/1000, 1/250 and 1/125 from right to left. The exponent r is given 
by the negative slope. 
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FIG. 6. Cluster radius R{s) in 2 to 6 dimensions (from left to right). The inverse slope yields 
the fractal dimension jj,. 
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FIG. 7. The correlation length ^ as function of (f/p) ^- The slope yields the critical exponent 
v. (□ = square lattice, A = triangular lattice, * = next-nearest-neighbor interaction.) 
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FIG. 8. Crossover scaling function F(x) for the correlation length for different values of the 
immunity. The dashed line represents i^(0) as obtained at ^ = gc- 

33 




FIG. 9. Dynamics on a string of A: = 4 sites. Trees are black, empty sites are white 
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FIG. 10. Size distribution of the fires for d = 1, //p = 1/25000 and L = 22°. The smooth hne 
is the theoretical result which is valid for cluster sizes < Smax- 




S 



FIG. 11. Normalized size distribution of fires in a version where the life time distribution of 
empty sites is a constant function P{t) = 1/100 for times r < 100. The parameters are L = 10000 
and / = 0.00001. The smooth line is the theoretical result, which is valid for s < Smax- 
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FIG. 12. Normalized size distribution of fires in the model with tree conservation, for p = 0.43 
and L = 512, 1024, 2048, 4096, (a) before and (b) after rescaling. 




FIG. 13. Stationary state in the model with tree conservation for p = 0.50, L = 2048 and 
absorbing boundary conditions (trees are black, empty sites are white). 
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TABLES 

TABLE I. Numerical results for the critical exponents in 1 to 8 dimensions (* = with logarith- 
mic corrections, ^ = calculated from scaling relations), taken from [26]. The exponents with index 



"perc" are those of percolation theory [33] . 



d 


1 


2 


3 


4 


5 


6 


7 


8 


L 


220 


16384 


448 


80 


32 


20 


12 


8 


T 


2 


2.14(3) 


2.23(3) 


2.36(3) 


2.45(3) 


2.50(3) 


2.50(3) 


2.50(3) 


Tperc 


2 


2.05 


2.18 


2.31 


2.41 


2.5 


2.5 


2.5 


X 


1* 


1.15(3) 


1.30(6) 


1.56(8)t 


1.82(10)t 


2.01(12)+ 


2.01(12) + 


2.01(12)+ 


1/5 


0* 


0.48(2) 


0.55(12) 












Pi 


1 


0.4081(7) 


0.2190(6) 


0.146(1) 


0.111(1) 


0.090(1) 


0.076(1) 


0.066(1) 


M 


1 


1.96(1) 


2.51(3) 


3.0 


3.2(2) 








Mperc 


1 


1.90 


2.53 


3.06 


3.54 


4 


4 


4 


U 


1* 


0.58 


0.52(3)+ 


0.53(3)t 


0.57(7)+ 








v' 


1* 


0.58 


0.64(6)t 


0.78(8)t 


0.92(10)+ 


1.04(11) + 


1.05(11)+ 


1.05(11)+ 
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